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Abstract. Giving new insight for line broadening theory for atoms with more structure 
than hydrogen in most stars. Using symbolic software to build precise wave functions cor¬ 
rected for 5s, 5p quantum defects. The profiles obtained with that approach, have peculiar 
trends, narrower than hydrogen, all quantum defects used are taken from atomic database 
topbase. Illustration of stronger effects of ions and electrons on the alkaline profiles, than 
neutral-neutral collision mechanism. 
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1. Introduction 

Here we define what is needed for the theory: We introduce the profile function normal¬ 
ized: F{A(o) requiring normalization: 

(1) / F{Aco)d{Aco) = 1 

J —oo 

The radiant power function, that is the power emitted thoug unit frequency is: 

4®^ 

(2) P{A(o) = -^F{A(o) 

We recall the two relations defining the Fourier transforms variables : (O and s (time vari¬ 
able) such that (O.s is dimensionless. 

/ oo 

F^“-^F{Aco)d{Aco) 

(3) /(A®) = - r <P{s)e‘^^-^ds 

n J-oo 

2. Hydrogen facts 

Recalling some facts on Hydrogen lines broadening: such as Balmer lines Ha, and 
Hj3 or Lyman lines : Lya, and Lyj3 hcoo being the energy gap between levels defining 
a transition i —> / the detuning A® giving rise to F{A(o) for hydrogen lines is mainly 
first order Stark effect: for instance exists for Hj3 whose wavelength at the centre of the 
= 4 —>■ n,' = 2 transition is: Xq — 486.Inm. 

These Hydrogen lines are broad because of the degeneracy of levels defined by \nlm > 
there can exists g = 2 x (2Z -f 1) for an 1 defined state such as 0 > Z < n — 1 for an identified 
H line such as Hj3 = 4 —> n, = 2 line we have N = 4.n}xnj that is 512 sublevels implied 
and 128 neglecting the electron spin. Good literature exists for such matters, that is taking 
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into account for peiturbers (ions and electron) effects on the radiating atom. [Feautrierl 
( [T976] |. 


(4) 

A£ = 

3 T 

-n X («! — n 2 )eF + 0{F^) 

(5) 

AE = 

hAco = h{(Oo - (o) 




(6) 

AE = 



These two approaches gives rise to the complete theory of broadening: 

1) when the static held 1^1 = Cst is constant relatively to decay rate, in fact a held generated 
by static heavy charges ions or protons (relatively to the electron)). Then the Stark prohle 
has such dependence: 


(7) /(A):^|AA|^ 

2) The held F(t) varies quickly before the atom relaxes, 


(8) 

r2(f) 

(9) 

r(f) ~ h + v.f 


The quantum theory of one electron +H (or radiating atom) system fVan Regemorter| ( [T972j l 
implies the relation tiAo) = ^ (k? — k^) is the good answer to take into account for large 
Ao) and replaces the impact parameter semi-classical approach. 


3. Alkaline lines 


Dealing with atoms having a structure such as alkalines: Li, Ca, K, Mg, Na or atoms 
such as He, and O. These atoms can be modelled using what is called since a long time [] 


Born (19261 the quantum defect. The optical electron : the one that gives rise to transition 


(quantum jumps (a —> j3) suffers from an additional potential: the polarization potential: 

CtD being the dipolar static polarizibilty. The full theory begins with recent 
review work |Schwerdtfeger| ^006) 1 and displays the development of an energy level as 
orders of the held strength F : (if F{t) = Eq) this the well known Stark effect. 


EiF) =Eo + ^dF + i ^dF^ + 0{F^) 

The polarization potential varying as is easily reckognized as the second order term 
dF^. Let’ us introduce the way to deal with the Vp{r) potential. 


(10) 

P = 

Ud-E 

(11) 

II 

—p.dE 

(12) 

II 

—aoE.dE 

(13) 

Vp{r) = 

ao 

2r4 

(14) 
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3.1. Semi-classical expression. How to deal with; Let’us start with the semi-classical 
formula for the prohle F{co) 


(15) 


Fi(o) = 


lim 




2kT 


This equation is put forward in Van Regemorter ( \912\ in his review of spectral line broad¬ 
ening. I adapt that question to a quite similar way: I need not have |'T;(f) > as a time 
dependant wave function, but the modifyed radial |'T(r) > with no time dependence. Giv¬ 
ing for alkaline species, with known quantum defects (there are data from Topbase): 


(16) 

(17) 

(18) 
(19) 


(X — 

p = nfjf^ 

\'Yp{ 7 )>= R^xYp{e,(t,) 


or for Hydrogen the well known basic ket |R„/(r) x > 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


yil - 

\^„i{r)>= RnixYU9,p) 

n'l' = Uflf 

>= ^nl' 


Our purpose is to dehne the most efficient way, the |'Ta(r) > using true quantum defects 
leading to n* ; the effective quantum number. Each atomic species has its peculiar quan¬ 
tum defect. These are now available from data base such as Topbase. Now it is a fact that 
the 1 kinetic momentum is dehned by eigen value of the spherical harmonics for a pure 
Coulomb potential, becomes /* = / — 5^ the degeneracy of the levels disappears. The quan¬ 
tum number set is :a = n* = n — 5/, =1 — 5/. The physical effect produced can be explain 

this way: the optical electon getting away from the closed shell beneath polarizes the core 
shell, the higher the levels of the optical electron , the nearest to hydrogenic ’’states” are 
the transitions. 


(24) 


L= 0 = S^ 5s 

L= l=P^5p 

5/ = 5i > 5p > 5rf > 5/ —> 0 


4. Time dependent method 

It is seen that one can change the semi-classical formula, into the way suggested by 
Schiff|(|1968|l, with no less generality. 

'Em 

'Tc((gO = 'Pa(GO) X e^ir- the same for the |j3 > | that is: 'i‘p{r,t) = 'T| 3 (r,0) x . 
The transition probability is proportional to ; 

I < > P X From the text book of L.I. Schiff, we use his g(t) plateau 

function to transform and the relation T = 1—5; 

(25) <p\iS-l)\a>=-j^<a\T\p > 

The P{A(o) function is then proportional to; 

Wap°^\< CC\T\P > p the transition probalibity . 
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The following operators are the same: 
A{t) = I < ^a{r,t)\D\^p{r,t) > ^ and 
B{t) = \<^'a{r)\D\'i>p{r>)\^xgit,to) 



ALKALINE LINES BROADENING IN STARS. 


5 


5. Toward the Mg profiles 


It is clear that one can obtain the wave functions of these alkaline elements such as Mgl 
neutral with some quantum defects; 

5, = 1.52 and 5p = 1.04 and 5d = 0.56 S=0 (Singlet) 

1 (Triplet) see Ref. 


5, = 1.63 and 5p = 1.12 and 5d = 0.17 S 


Born 


(1926l,ppl90 


There the Ik is the Kostelecky index , taking values such as 1 or 2 jKostelecky (19841 


This is done to insure the positivity of quantum numbers nj^. For a Mg3s 
Mg3s —> 4p(5 = 0) , the quantum numbers are below: 


3p or a 


(26) 

a = 

00 

II 

1 

CO 

II 

.* 

(27) 

//'* = 

0 — -f Ik 

(28) 

l'I'«(p) >= 

^a(p) X 0) 

(29) 


rift = 3 — dp = 1.96 

(30) 

h* = 

l-5p 

(31) 

|'F^(r)>= 

Rp{r)xYp{0,(p) 


The full theory of the calculation of wave functions of such type is done in |Bates| ( |1949| l, 
it implies subtle transformations with the theory of the Kostelecky index, related to super- 
symmetry transformations Kostelecky ( |1984| l. The quantum analog to the classical Born 
atomic model, with the precession of the ellispse, is obtained by |Bates| ( |1949[ l. The good 
quantum theory needs to consider a modification of the T/m(0,0) —> 0) It is useful 

to compare two transitions such as ; Mg3 As triplet quantum defects to consider whose 
wavelength is A,/ = 517.83nm with an Hydrogen line such as Hhp 4s A,y = 1875.1 l«m. 
All data taken from NBS (19691.The ratio of the intensities is given by; 

I < anH.,(H.k.cos(0)|j3nl,l' 


(32) 


R = 


> 


'Lm=-l I < n,l,m, |r.cos(0)|«',/± l,m > p 


The ratio R can be considered as the ratio of the Einstein An- of the two lines, for such two 
lines ; R = 0.56 x 10^® 

For most ionic species, whose single or optical electron gains high distances 
< r >— rv^QQ > 20ao, the quantum defects disappear, leading to simple hydrogenic be¬ 
haviour, the wave functions turn to be these of Hydrogen. 


5.1. Toward the global theory. We use the upward definition for the line shape; 


(33) 

F{(0) = 

J —oo 

(34) 

8it) = 

m-m 

(35) 

F{co) = 

\<^a\D\'Vii>\^ J\-‘^'.g{t)dt 

(36) 




These Fourier transforms are easily performed with ’’Mathematica” , and I write here the 
g{t,to) Jo being the width of the ’’plateau” in seconds. 

(37) g{t,to) = — fo) -!-//(—f -ffo) + H {t — to) + H {—t +to) 

Here is the picture of the distribution function g{t,to): 

It is reallly interesting to have fast symbolic software such as ’’Mathematica”, to perform 
the Fourier Transform of the function g{t,to) assuming that fo > 0 to get F{co,to), and have 
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negative second way to g(t) Schiff funetion 



Figure 1 . g{t,to) function to = 3s fo = 5s 


it as a good analytical function: 


(38) 


F{co,to) = 


/ia/ 2(1—e 'o®+ 2®cos(®fo)) 
i/tt® 1 + ®2 


It is interesting to look at the behaviour of the profile function F{(0,to) , this function 
oscillates because of the phases contained in the expresssion. 

For a given ® the shorter is the to width the least the profile F{(o) oscillates. To obtain 
the good former profile function we need to replace in the defined function F{(0,to) the 
variable ® by the variable A® = ® — ®o. Here are some pictures of the alkaline Mg 
element: 

Here ®o is the frequency associated with the center of the line ®o = —. 
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The parameter /i = | < '?«> |^ . should be the intensity at the center of the line that 
is 0) = 0)0 .That definition of Ii shall be modified to take into account the pressure effect. 


6. The cut-off radius and its significance upon the profile /^(AoTo). 

A simple and efficient way to take into account the pressure effect whose origin comes 

from external parameter such as density Na and thus the cut-off is the following La = . 

This cut-off parameter La is the one to use when there are no charges, and the surrounding 
of the emitters is mainly done of atomic species of density Na- If the emitters of the 
transition a —?► are part of a plasma, with a degree of ionization: T = 
the good cut-off could be the shorter of the the two Stehle ( 1996|l that is : 


(39) 

XDe = 

1 kTe 

V 4ne^.Ne 

(40) 

XDe = 

X 

o 

oo 

(41) 

XDi = 

XD, 


Pc = 

Min{XD{Ne,T),LA) 


XDN{ions) that is the Debye radius to compare the spatial extension La, because of the 
distribution of the charges in the plasma. 


7. Effect of the cut-off on the value obtained for the dipolar squared 

MATRIX ELEMENT ^fa/3 = l < 0:|D|j3 > ^ 

That simple explanation to take into account for pressure effects on radiating atoms 
,parts of such media, such as exists in Astrophysics: -Star atmospheres, -Dust in interstellar 
matter. -Molecular clouds. Finding the way to define the upper limit pc for the emitting 
atoms ,and to build the wave functions with their peculiar quantum defects for each species 
(Li, Na, O ,A1 , Ca ,Na ,Mg and even He) now easily reached with modern symbolic 
calculation software, I define there the a probability function p (v): 


(42) 


PapixJ) 


\jfUa\r\l5>\^ 


8. Some RESULTS coming FROM THE Pap (Xj T) function 

The variation domain of this function is the-half real axis , strictly positive as suitable 
for a length. 

We can consider such function: 

(43) 0<Pa.p{x{T)<l 

Pap (0) = 0 void of particle. 

Pap (-t) = 0 squeezed state of the optical electron. 

Pap (°°) = 1 existence of the particule. 
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9. Analytical results for F{A(o,to) taken from 

In such matter we delt with; two distincts problems are solved; The non hydrogenic be¬ 
haviour of the emitters or absorbers (He, Li,0 ,Mg ,Ca ,Na ,K) and their ionized species, is 
modelled with the quantum defects from which exist a large litterature and even databases 
such as Topbase. The correct building of the wave functions is a solved problem. |de Ker-| 
Itanguyl l flQTQ | l. The second problem solved in this paper, is the good way to obtain the 
profile F{A(o) from the Fourier transform •!>(?), whose value includes the restriction of the 
wave function through the cut-off parameter, pc = x{T). As a matter of fact an oscillation 
effect under the envelop of the F{A(o) < 

being the limiting impact profile function for light perturbers as electrons. There the 
transition frequency oio for a ns{5s) mp{5p) transition is given by; Ih = 13.606ey and 


e=\q\ = 1.6021910- 

*®Cand AEpip 

is given in atomic units ua. 

(44) 

®0 = 

2.481810''^7rx/// X A£„^ 

(45) 

®o = 

e X (Oq{Hz) 

(46) 


--) 

2 n% m% 


n* = 

n-5s 


OT* = 

2 

1 


Setting A® — co — coq. 

Fixing the fo parameter that is the width of the smooth function g{t,to). 

I rewrite down the Fourier Transform with the modified constant Ii changed into a F;^/- 
resulting from the shortening of the radial integration domain or ’’incomplete” dipolar 
squared matrix element. 


rX 

OpeQDFunction(x,ind,3,4) = / (f)^n= 4 -S (f)dr 

Jo ■' ’’ 

ao is the Bohr radius. 

(47) Ii=nlx {OpeQDFunction{0.l,l,3,4)^)aQ 

The Ti^f parameter is allways defined by the chosen transition / —> / is obtained by solv¬ 
ing; 

(48) g{t,to) X = OpeQDFunction(t , l,n,',ny-) 


(49) 


F(A®,fo) = 


_ r,_,y's/2 (1—e ^Oco+2c!)cos{w/o)) 


v/tCO ^ 1+0)2 

This function F{Aco,to) can be integrated on the A(0 variable from > oo. Using the 


Stehle 


(1996 1 , as a reference a function such as FI (A®) = 


(Ao))^ 


can not be normalized 


when integrated from —oo —oo while a function such as F(A®,fo) can. 
In fact; 


(50) 

(51) 


Fl{{Aco)dAco = 
F{{Aco,to)dAco = 


0 
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10. Conclusions 

This work whose conclusion is a new insight to complex atoms (more structure than Hy¬ 
drogen) line broadening, lies on two different points. First, the pressure effect is included 
in the calculation of the probability function Pap{x{T,N) , where x(N,T) is the cut-off pa¬ 
rameter Pc to be used. The surrounding of the radiating atom will fix it. It is clear that 
when x{N,T) —> the isolated atom line profile is obtained, in between the pressure effect 
is taken into account!. The second interesting fact is included in the definition of the width 
to of the smoothed g{t,to) distribution function: It can exhibit oscillations of the profile 
F(Aco,to). The choice of the width to should be driven by a sound consideration of the 
surrounding atoms or ions background. It seems plausible the densest is the medium of 
the surrounding plasma or bath of neutral perturbers, the shorter the to should be! Here are 
some sketches in Fig.4 and Fig.5 of the F{A( 0 ,to) for different width to whose range goes 
from to = Ips to to = 20i 
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Profile function F(Anj,tO) Il=|<Q'|D]yS>p=1 


2. X 10“^^ 


-2.x 10“^^ 
-4. X 10“^^ 


-5 X 10 


5x10 


-1x10*’ 


-1.03091 X 10*' -1.03091 x 10 


1x10* 


Figure 2. 3D Plot ofF{AcoJo) = profile func¬ 

tion (Oq = 1.0310*^. 


It happens that if Na — Ne the cut-off radius is allways smaller ta = .This is 
a proof that the pressure effects in a plasma with a charge density Ne is bigger than the 
pressure effects due to atoms of same densities Na = Ne- 
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Na 

IxlO^* 


4x10^' 6x10*' 

Long range radii 


O' 

2x10-' 4x10-' 6x10-' 8x10-' 1x10" 

Figure 3. dark blue line: cut-off radius in Bohr unit as func¬ 

tion of the atoms density Na in cmT^. light blue curve : cut-off radius 
Tg = 6.9 X 10® for an electronic temperature of Tg=\000K and a 
density of charges Ne- 
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Figure 4. Oscillations of the F(Aa),fo) function for fo l/r , 10/rs, 10 ms 
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t3=l CosTF 





Figure 5. Oscillations of the F{Acojo) function for to , 0.01 ms, 0.1 s, 
Is, 7 s ,20 s and all profiles 
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5 

4 

3 


' .... 

6x10^' 8x10^* lx 10” 

Figure 6. Display of cut-off bounds in unit of principal quantum num- 

_ 1 

ber two approaches L = yellow compared to XD Debye shielding 
radius brown. Explanation ; for a density of N = Na = Ne = 410^*cm^^ 
the low n=2 is obtained for Xd and n=4 for La — 3 .The cut-off pc 

varies from 4ao to 16ao. from charged to neutral surroundings. 

T=100,500,100,5000,20000 K 
5 - 

4 - 

3 - 


2x10-' 4xl0-‘ 



2xl0^‘ 4x10^' 6x10^^ 8xl0^‘ 1x10^^ 


Figure 7. Debye length XD{Na, T) for several temperatures 













